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In this work we reproduce tiie piienomenology of tiie electromagnetically induced transparency 
and dynamic Stark efltect in a dissipative system composed by two coupled bosonic fields under linear 
and nonlinear amplification process. Such a system can be used as a quantum switch in networks 
of oscillators. 
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Much of the recent work in quantum optics has been 
devoted to quantum coherence and interference as a re- 
sult of the mastery gained in manipulating atom-field 
interaction, which allows us to probe fundamental quan- 
tum effects. Within this scenario, electromagnetically- 
induced transparency (EIT) plays an essential role 
in a variety of processes, ranging from lasing without 
inversion 2] and enhanced nonlinear optics [Sj to quan- 
tum computation and communication Relying on 
destructive quantum interference, EIT is a phenomenon 
in which the absorption of a probe-laser field resonant 
with an atomic transition is reduced or even eliminated 
by the application of a strong driving laser to an adja- 
cent transition. This effect can be understood as arising 
from the ac-Stark splitting of an excited atomic state 
linking the two adjacent transitions. When the splitting 
is smaller than the excited state width, the two result- 
ing levels are indistinguishable, leading to a destructive 
interference in the probe absorption spectrum 0. As 
the intensity of the driving field increases, the splitting is 
enlarged, eliminating the indistinguishability and, conse- 
quently, the absorption spectrum evolves to the Autler- 
Townes (AT) doublet. By a similar mechanism, although 
relying on nonlinear interactions, the strong excitation of 
a two-level atom induces a dynamic Stark effect (DSE) 
leading to the appearance of sidebands in the resonance 
fluorescence spectrum @. 

Recently, some effects of three-level atom optics, in- 
cluding EIT, were simulated |6| via the tunneling interac- 
tion in the context of trapped atoms. In this connection, 
the phenomenology observed in EIT was reproduced in 
a system of two coupled dissipative classical oscillators 
in Ref. Q. One of these oscillators, modeling the three- 
level atom, is subject to a linear driving force 7] playing 
the role of the probe field. The pump field is simulated 
by the coupling between the oscillators. The quantum 
version of this classical analog of EIT, which helps to 
deepen our understanding of this phenomenon and its 
properties, will be analyzed in the present work. Here, 
the classical oscillators are replaced by two quantum res- 
onators interacting through a Josephson-type coupling. 
Linear and nonlinear amplifications are employed as the 
probe fields, the former is used to emulated the EIT phe- 



nomenon, and the latter is applied to obtain the DSE — 
reproducing the phenomenology of the nonlinear mecha- 
nism related to resonance fiuorescence from a strongly- 
driven two-level transition. The coupling parameter be- 
tween the resonators plays the role of the driving field, 
by analogy with the standard EIT and DSE. We show 
that the signatures of these processes can be obtained by 
measuring the field quadratures. 

The Hamiltonian of the coupled dissipative resonators, 
labeled hy t = 1, 2, is given by {h = 1) 
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where a\ (ag) is the creation (annihilation) operator for 
the oscillator mode i of frequency u!(, whereas fojij, {b(k) is 
the analogous operator for the kth bath mode of oscillator 
£, whose corresponding frequency and coupling strength 
are ujik and Vik, respectively. The coupling strength be- 
tween the oscillators is A and the classical driving field 
applied to oscillator 1 has intensity F and frequency pv, 
with the parameter p = 1 or 2 for the linear or nonlinear 
amplification process. 

Following the reasoning in Ref. P| , we obtain (in a ro- 
tating frame with frequency v and assuming the reservoir 
at zero temperature) the master equation for the system 
dynamics, given by 
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where is the damping rate of resonator i and we have 
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defined 
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with Qi = LUg ~ i> being the effective frequency of the 
oscillator £ shifted by the amplification process. Using 
the standard procedures and assuming Markovian white 
noise, we derive a c-number version of the master equa- 
tion ^ for the symmetric ordered characteristic function 
X{{ve},t), given by 



dt 



E 



+ir(r7i+??i*)(»7i-'7i*)*'' ^Uti 



(4) 



where C^({?7f}) = + + 2« (p - 1) F -ql, 

= + iil^, and is the complex parameter of 

the characteristic function. From the analytical solution 
of Eq. , obtained by standard techniques @ , we ob- 
serve for the linear amplification (p = 1), that the modes 
in both resonators always evolve to stationary coherent 
states, characterized by the equilibrium between the am- 
plification process and the dissipative mechanisms, medi- 
ated by the coupling strength A. These stationary coher- 
ent states, \a)-^ and |/3)2, have the complex amplitudes 
a = -iFk^l (A2 + A1A2) and /3 = -\F I {>? + A1A2). 

The EIT-like phenomenon, which occurs in the re- 
gion where Fi ;» A ^ is characterized by the 
state of the field in resonator 1 around the relaxation 
time of the joint system. For the linear amplifica- 
tion case, the relaxation time, derived from the condi- 
tion 9x({?7f} ,t)/dt 0, is given by the simple ex- 
pression Tfl = 2AVri + r2/2. In Ref. @], it is noted 
that the cavity-field states are interchanged between the 
resonators through state-swap and the state-recurrence 
dynamics. Before the relaxation time, the coherence dy- 
namics of quantum states between the two resonators 
prevents the EIT from occurring. Therefore, the system 
must achieve a stationary state to enable the appearance 
of the destructive interference, which promotes the EIT. 
For this reason, from here on we will be interested in the 
states of resonator 1 around the relaxation time. 

In atomic samples the EIT phenomenon is character- 
ized by the permittivity of the medium, associated with 
the complex quantity n -f ik, where the dispersive and 
the absorptive responses follow, respectively, from the 
real part of the refractive index n and the absorption 
coefficient k 0. Here, the role of the permittivity is 
played by the quantity S^^'^ e*®* representing the center 
of the quasi-probability distribution of the cavity-field 
state in the phase space defined by its quadratures Xi = 



+ aj^ /2 and Yi — — aj^ /2i. The absorption 
coefficient is directly associated with the mean energy 
of the cavity field given by the correlation function 
£g = l^Xj + (in units of huj(). The dispersive re- 
sponse is associated with the angle <d through the corre- 
lation function {XgYt -f YiXi) jZi = sin (20^). 

In Ref. the dispersive response Re[a;i(cj — Si)] is 
derived from the real part of the frequency dependence of 
the amplitude x\i^ — of one oscillator (frequency wi), 
submitted to a linear driving force (frequency w), and 
coupled to another oscillator (frequency Zj^ and ampli- 
tude X'l). In the present model, the equivalent mean value 
(l^) also leads to the required dispersive response, as in 
[3, for the linear amplification process (p = 1). However, 
in the case of nonlinear amplification (p = 2), the mean 
value (y^) is not a convenient quantity to measure the dis- 
persive response since it is null when the initial excitation 
of the cavity field I is also null (vacuum state). To cir- 
cumvent this problem, we have considered the correlation 
function (XtXi -|- YgXi) /8i to describe the dispersive re- 
sponse in a quite general way. In fact, its classical analog, 
proportional to Re [x^pg] (where pi is the momentum of 
the classical oscillator £), also leads to the usual disper- 
sion exhibiting a very steep slope in the EIT regime. The 
absorptive and dispersive responses given by the corre- 
lation functions -|- F/) and {XgYg A-YgXi) / £1, re- 
spectively, are computed from the characteristic function 
x{{Vm}.t) (to = 1,2) 
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Linear amplification. To characterize the EIT, we plot 
the mean values £1 and sin (20i), for resonator 1, against 
the detuning A = {v — loi) /Fi. Before the relaxation 
time {tr), i.e., t <^ tr, the EIT is inhibited by the co- 
herence dynamics between the cavity-field states. Near 
the relaxation time, the EIT starts appearing, evolving 
to its standard shape at the relaxation time, as displayed 
in Fig. 1(a). In Fig. 1(b) we assume the non-degenerate 
case LUi = LU2 — £ But the relaxation time, where the hole 
burned when wi in the degenerate case [Fig. 1(a)], 
shifts to ~ LUi + e (where e/Fi = 3). It is interesting 
to note, that in the non-degenerate case [Fig. 1(b)], the 
EIT happens when the linear amplification (h') is in res- 
onance with resonator 2 (^2). A complete transparency 
occurs only when r2 = 0, since the minimum of the hole 
burned in the peak of the absorptive response is given 
by ff*" = [2rr2/ (4A2 -hrirz)]^ Assuming the same 
parameters as those in Fig. 1(a) but a stronger coupling 
strength (A/Fi = 1), the EIT evolves to the AT regime 
as shown in Fig. 1(c). Finally, the dispersion curves, in 
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FIG. 1; Absorptive-like (left side) and dispersive-like (right 
side) curves for the linear amplification process (at t = tr) 
assuming: (a) the degenerate case {ui = UJ2) and r2 <C A <C 
Fi, (b) the non-degenerate case (uj\ = uj2 — e),. (c) the de- 
generate case and F2 ^ A « Fi. To plot these figures we 
have assumed the following parameters: F2/F1 — 10"'' and 
r/Fi = 1/4. For (a) and (b) we have used A/Fi = 2 x 10"^ 
and for (c) A/Fi = 1. 



the right side of Fig. l(a,b,c), display the weU-known 
signature of the EIT effect. 

Nonlinear amplification. In the nonhnear or paramet- 
ric amplification case (p = 2), there is a well-known 
threshold in the dynamics of the system 

Ea, JJJ which 

is reached when the rescaled amplification parameter, 
^ = 4F/ri, equals the critical value (which depends on 
the critical amplification amplitude F c)- We obtain three 
different dynamics for the solutions of Eq. which 
depend on whether the rescaled amplification is strong 
> ?c), critical (^ = ^c), or weak (^ < ^c)- When 
^ > the asymptotic solution diverges for any physi- 
cal quantity whereas we obtain stationary solutions when 
^ < ^c- In the degenerate resonance case, loi = UJ2 = 
the critical amplification parameter is given either by 

= 1 + 4AV (FiFa) (if A < T2/2) or = 1 + r2/ri(if 
A > T2/2). Evidently, when A = we recover the usual 
result = 1 for an uncoupled system under nonlinear 
amplification. In these three regimes the state of oscilla- 
tor 1 evolves to a squeezed coherent state (SCS) which is 
stationary for ^ < and whose energies increase mono- 
tonically for ^ > ^c- 

Now, in order to analyze the behavior of the absorp- 
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FIG. 2: Absorptive-like curves for the nonlinear amplification 
process {ui = LJ2): in (a) the weak regime (^ = ^c — 1.5x 10~*, 
= 1 + 10"'' ), (b) critical regime (^ = ^c), (c) strong regime 
(5 = + 1,5 X 10~*), and (d) in the weak amplification 
regime at the relaxation time tr, but increasing the coupling 
strength. To plot these figures we have assumed the following 
parameter: r2/Fi = 10~*. For (a-c) we have used A/Fi = 
5 X 10"^ and for (d) we have used A/Fi = 10"^ Where the 
dotted line is for t — 0.8 x tr, the solid line is for t = tr, 
the dashed line is for t = 1.1 x tr, and we have assumed the 
vacuum states as initial states for the resonators 1 and 2. 



five and dispersive curves of resonator 1, in the three 
regimes exposed above, we plot the rate £i/£, instead of 
£1 , where £ stands for the value of £1 when the detuning 
A = 0. In this case, the EIT occurs only in a narrow in- 
terval 1^1 ^ + 10"* around the threshold given by the 
critical parameter ^c, i.e., aroimd Fc — (Ei + 4A^/r2)/4 
for A < T2/2, or Fc = (Ei + E2)/4 for A > E2/2. Dif- 
ferently from the linear amplification process, where the 
EIT can be achieved independently of the ratio F /Fi, 
the strong dependence of this effect upon the amplitude 
of the nonlinear driving field follows from the two-photon 
nature of this amplification process. As far as the cou- 
pling strength is concerned, although it must be much 
smaller than Fi as in the linear amplification case, here 
it can be around E2, obeying the relation E2 < A ^ Fi. 

In the nonlinear case the relaxation time th can only 
be obtained from the condition dxUve} 1 1) / dt — > in 
the weak amplification regime. This time will be used 
as a reference to compare the three different amplifica- 
tion regimes, since there is no stationary dynamics when 
C > Cc- We can see in Fig. 2(a), the absorptive curve for 
the weak amplification regime, that the EIT is not so pro- 
nounced as in the linear amplification case. In fact, it can 
barely be characterized as an EIT, since the maximum 
of the ratio £i/£o is close to unity (i.e., the minimum 
of the hole burned in the peak of the absorptive curve, 
when A = 0). The absorptive curve remains immovable 
for t > Tji due to the stationary solution achieved in the 
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weak amplification regime. Differently from Fig. 2(a), we 
observe that in the critical regime, Fig. 2(b), the maxi- 
mum of the ratio Ei/£ is about two orders of magnitude 
larger than unity {£i = £) for t = T]i (solid line), increas- 
ing as time goes on. However, as in the weak amplifica- 
tion regime, the behavior observed here can hardly be 
characterized as an FIT, but rather a mixing of FIT and 
AT effects, since the hole burned in the peak of the ab- 
sorptive response presents a significant width compared 
to the usual FIT. In Fig. 2(c), for the strong amplifica- 
tion regime, a mixing of the FIT and the AT effects is 
observed again, but the ratio £i/£ does not increase as 
quickly as in the critical regime — even so, £i and £ in- 
crease, independently, as faster as they do in the critical 
regime. The maximum of the rate £i/£ is (only) about 
one order of magnitude larger than unity {£i — £) for 
t = tr and increases as the time goes on. We note that 
the range of the detuning A, in Figs. 2(b) and (c) — 
where the absorptive curves assume significant values — 
is considerably smaller than that in Fig. 2(a) for the 
weak regime. Therefore, as the rescaled amplification pa- 
rameter ^ increases, the range A of the absorptive curve 
decreases. From this perspective, it is not difficult to un- 
derstand why the range A of the absorptive curve in the 
linear amplification regime (Fig. 1) is larger than those in 
the nonlinear one. Regarding the dispersive curve, it hap- 
pens to be practically the same for the three parametric 
amplification regimes and the steep slope characteristic 
of the FIT effect is lost. 

Finally, in Fig. 2(d) we display the absorptive curve in 
the weak amplification regime at the relaxation time tr, 
but increasing the coupling strength such that A/Fi = 
10~^ (instead of 5 x lO^'"^). We thus observe the well- 
known pattern exhibited by the resonance fluorescence 
spectrum of a two-level atom driven by an incident field 
whose Rabi frequency is comparable to, or larger than, 
the atomic linewidth. After all, as the phenomenology 
of the FIT can be reproduced in a system of two cou- 
pled dissipative oscillators, the same occurs with the phe- 
nomenology of resonance fluorescence. Resonator 1 again 
plays the role of the two-level atom, while the strong 
driving field (whose intensity leads to a modulation of 
the quantum dipole moment inducing sidebands in the 
spectrum |12)|) is simulated by the nonlinear amplifica- 
tion together with resonator 2. In fact, high- intensity 
driving fields give rise to nonlinear interactions which 
may be attributed to two laser fields modeled by 
the actual parametric field plus resonator 2. The pat- 
tern displayed in Fig. 2(d) can also be obtained in the 
critical and strong amplification regimes but, of course, 
it does not occur with linear amplification, whatever its 
intensity. The parametric amplification field is the only 
nonlinear ingredient in our system of two coupled dissi- 
pative oscillators. 

We have shown that a system of two coupled dissi- 
pative oscillators, under linear and parametric amplifi- 



cation, can reproduce the phenomenology of FIT and 
DSF, respectively. Such phenomena can be, at principle, 
implemented in several physical contexts, such as: cav- 
ity quantum electrodvnamics |l4l| , trapped ions [isj , and 
nanomechanical oscillators [ig]. The system presented 
here can be used for quantum switching purposes, since 
the hole burned in the FIT regime can be shifted by ex- 
ternal means, feeding or draining the field in resonator 1. 
For example, the dispersive interaction of a fermionic par- 
ticle with the field in oscillator 2 can pull both oscillators 
out of resonance, changing the profile of the absorption 
curve of oscillator 1 from Fig. 1(a) to 1(b). Therefore, 
these effects have the potential to be used as a control 
mechanism in networks of quantum oscillators, a subject 
which attracts great attention nowadays [l^ . 
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